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Abstract 



The ratio of the CKM quark- mixing matrix elements |V^,|/|Vrf,| has been mea- 
sured using B hadron semileptonic decays. The analysis uses the reconstructed 
mass Mx of the secondary hadronic system produced in association with an 
identified lepton. Since B — > X u lv transitions are characterised by hadronic 
masses below those of the D mesons produced in B — > X c tv transitions, 
events with a reconstructed value of Mx significantly below the D mass are 
selected. Further signal enrichments are obtained using the topology of re- 
constructed decays and hadron identification. A fit to the numbers of de- 
cays in the b — > u enriched and depleted samples with Mx above and be- 
low 1.6 GeV/c 2 and to the shapes of the lepton energy distribution in the B 
rest frame gives IK&I/IK&I = 0.103±g;g^ (stat.) ±0.016 (syst.) ± 0.010 (model) 
and, correspondingly, a charmless semileptonic B decay branching fraction of 
BR(5 -> XJLv) = (1.57 ± 0.35 (stat.) ± 0.48 (syst.) ± 0.27 (model)) x 10~ 3 . 



(Phys. Lett. B478(2000)14) 
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1 Introduction 

The measurement of the branching ratio for the decay b — > uiv provides the most 
precise way to determine the \y u \\ element of the Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix. Evidence for a non-zero value of \V u b\ was first obtained JT],0] by observing 
leptons produced in B decays with momentum above the kinematic limit for b — > civ 
transitions. However, extracting \V u b\ from the yield of leptons above the b — *> civ endpoint 
is subject to a large model dependence. More recently, exclusive B — > uiv and B — > piv 
decays were observed and their rates measured But the determination of \y u }\ from 
exclusive semileptonic decays also has a significant model dependence. 

The extraction of \V u b\ from the distribution of the invariant mass Mx of the hadronic 
system recoiling against the lepton pair in B — > X u iv transitions was proposed several 
years ago || and it has recently been the subject of new theoretical calculations |||7fl. 
There have been two other \V u b\ determinations at LEP based on inclusive analysis of 
semileptonic decays ||. The method used here starts from the observation that in most 
B — > X u iv decays the hadronic system recoiling against the iv has an invariant mass 
below the charm mass (see Fig. [I]). Because a much larger fraction of the total rate is 
involved, the model dependence when extracting \V u b\ from the decay rate to such states 
is much smaller than when using the decay rate to leptons above the b — > civ endpoint 
or that to exclusive final states [f|[7|]. 

This paper presents the first determination of iV^I/jV^I based mainly on candidate B 
semileptonic decays with reconstructed hadronic invariant masses below the D mass and 
enriched in b — > u transitions using the secondary vertex topology and identified kaons 
and protons. The shape of the lepton energy spectrum in the B rest frame is also used. 

Section 2 describes the event preselection, the particle identification, the reconstruction 
of the hadronic secondary system and of the B energy and direction, and the b — > u 
enrichment. Section 3 presents the final event sample, the extraction of | V^l/I V^|, the 
stability checks, and the evaluation of the systematic errors. Section 4 summarises. 



2 Data Analysis 

The analysis was performed using data collected by the Delphi detector at Lep at 
centre-of-mass energies around the Z° pole between 1993 and 1995, corresponding to 
2.8 xlO 6 Z° hadronic decay candidates. The Delphi detector was described in detail 
in U and its performance was reviewed in JTIJ. The backgrounds were estimated using 
samples of Z° hadronic decays generated with Jetset 7.3 |TTJ and passed through the full 



detector simulation. These simulated events corresponded to 4.9 times the dataQ and were 
evenly divided in order to describe the Delphi detector response in the different years 
of data taking. Background B —>■ X c iv decays were generated using, for the exclusive 
modes, form factors based on a relativistic quark model [O. 

Events containing signal B — > X u iv decays were simulated using a dedicated decay 
generator JTB| interfaced with Jetset and passed through the full Delphi detector sim- 



ulation. Hadronic final states were produced using a tuned version of the parton shower 
model. The values of the branching ratios for the exclusive B — > uiv and B — > piv decays 
were forced to those measured by Cleo HH. The probability of producing vector and 
axial vector resonances was tuned to agree with the measurements of their inclusive cross 
sections in Z° decays. 

1 The numbers of simulated events quoted hereafter in this paper have all been renormalised by this factor 4.9 so as to 
be directly comparable with the numbers of data events quoted 
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2.1 Event Preselection and Reconstruction 

Hadronic events were selected using the standard Delphi criteria []14|| . These yielded 
2789419 events in the combined 1993-95 data. The sample was enriched in Z° — > bb events 
by applying a 6-tag algorithm based on measurements of the track impact parameters. 
This algorithm computes the probability for all reconstructed particles to originate from 
the event primary vertex ||15|| . This probability was required to be smaller than 0.08 
corresponding to an efficiency of 85% for Z° — > bb and a rejection factor of about 7 for other 
hadronic decays. Events were divided into two hemispheres by a plane perpendicular to 
the event thrust axis. Jets were reconstructed using the Luclus clustering algorithm [|Il[ 
with a dj oin value of 6.0 GeV. Only events with two or three jets were used. For the 
two most energetic jets in each event, a secondary vertex reconstruction in the jet was 
performed using those charged particle tracks with significantly large impact parameters. 
This procedure allowed to inclusively reconstruct a jet secondary vertex in 60% of all the 
jets in events satisfying the 6-tagging criteria. 



2.2 Particle Identification 

For this analysis, hadronic 6-tagged events were required to contain one identified lep- 
ton (e or /i). Candidate leptons from semileptonic B decay were selected in the momentum 
interval 3.5 GeV/c < p < 25 GeV/c. 

Muons were identified by the hits associated in the muon chambers. The efficiency 
was estimated from simulation to be (83.0±2.0)%. The probability for a hadron to be 
misidentified as a muon was measured on data to be (0.68±0.03)%. 

Electron candidates were selected using a Neural Network based on the response of the 
HPC electromagnetic (e.m.) calorimeter and on the measured specific ionisation (dE/dx) 
in the TPC The efficiency of this selection criteria was measured with Compton events in 
data, yielding (70.0±2.0)% with a misidentification probability of (0.38±0.03)%. In order 
to reduce the background from b — > c — > I and c — > i transitions the lepton candidates were 
required to have p\ n > 0.5 GeV/c, where p\ n is the momentum transverse to the jet axis 
reconstructed including the lepton candidate. In order to ensure its accurate extrapolation 
to the production point, each lepton candidate was required to have at least one associated 
hit in the silicon Vertex Detector and a positive lifetime-signed impact parameter relative 
to the primary vertex. 

The identification of strange mesons, which are produced in the cascade b — > c — > s 
decay, was used to reduce the background from charmed b decays. Kaons and protons, 
with p > 2.5 GeV/c, were identified by the combination of the response of the Delphi 
Ring Imaging CHerenkov (RICH) detectors and the dE/dx in the TPC [|l6j. K® were 



reconstructed in their 7r + 7r decay mode ||10|| , requiring a total momentum larger than 
2.5 GeV/c. 

Candidate 7r° — > 77 decays were also tagged. At energies up to about 6 GeV, the 
two photons from the 7r° decay are separated enough to produce two distinct electromag- 
netic showers in the e.m. calorimeter. Photon pairs with a total energy below 9 GeV 
and with an invariant mass 0.045 GeV/c 2 < m 77 < 0.225 GeV/c 2 were accepted as tt° 
candidates. At larger energies, 7r°'s were discriminated from photons by the e.m. cluster 
shape reconstructed in the HPC calorimeter. 
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2.3 Reconstruction of the Secondary Hadronic System 

A total of 52952 hemispheres with an identified lepton were accepted in data and 52661 
for the simulated backgrounds. 

The secondary hadronic system was reconstructed using an inclusive method. First a 
likelihood variable was computed as the product of the likelihood ratios that the particle 
originated from the B decay or from the primary hadronisation for a set of discriminating 
variables. These were divided into two categories. The first consisted of the following 
six kinematical variables: transverse momentum p t , ratio of particle momentum and jet 
energy p/Ej et , particle mass, rapidity, rank in decreasing energy order, and increase of the 
invariant mass of the particles associated to the jet secondary vertex by the addition of 
this particle. These variables were computed for all charged particles with momenta above 
0.5 GeV/c and neutral particles with energy above 1.0 GeV. The second category grouped 
four topological variables computed only for charged particle tracks with p > 0.7 GeV/c 
and associated hits in the Vertex Detector. These variables were the track impact pa- 
rameters relative to the primary event vertex in the R — <p and z projections normalised 
to their errors, the y 2 contribution of the track to the jet inclusive secondary vertex, and 
the distance of the point of closest approach of the track to the jet axis normalised to the 
distance between the primary and secondary vertex. 

The hadronic secondary system was then reconstructed via an iterative algorithm that 
used the charged particles sorted by decreasing value of the above likelihood variable to 
define a secondary vertex. All charged particles with likelihood larger than 0.75 were 
tested for their compatibility with originating from a common secondary vertex. Those 
contributing less than 4.0 to the x 2 of the secondary vertex fit and giving a total invariant 
mass Mx smaller than 3.0 GeV/c 2 were accepted. After the charged hadronic system was 
defined, identified K® and tt° candidates with a likelihood larger than 0.65 were tested. 
At most two 7r°'s were accepted at each vertex, if the total invariant mass of the particles 
associated to the vertex did not exceed 3.0 GeV/c 2 . Reconstructed secondary hadronic 
systems with only neutral particles, with an absolute value of the charge above 1, or with 
a total energy smaller than 4 GeV were rejected. Secondary hadronic systems consisting 
of a single particle were accepted if the particle was consistent with forming a common 
vertex with the lepton, in this case the mass Mx was set to the 7r mass. The total 
and charged multiplicities of the secondary system for data and simulation are shown in 
Fig. 0. The corresponding energy distributions for the secondary hadronic system X, the 
XI system, the missing energy in the hemisphere, and the Xlv system are shown in Fig. 0. 
A secondary hadronic system was reconstructed in 75% of simulated B — > X c lv decays 
satisfying the event selection criteria described above, and in 69% of signal B — > X u £v 
decays. 37986 decays were reconstructed in the data and 37899 in the simulation. 

The total mass of the candidate B decay, Mxe&, was estimated from the invariant mass 
of the system formed by the secondary hadronic system, the lepton and the neutrino as 
discussed below. For B — > X u £u and fully reconstructed B — > X c tv decays, M X ep peaks 
at ~ 5.0 GeV/c 2 and has a resolution of 0.9 GeV/c 2 , while for partially reconstructed 
B — > X c £u decays Mxtp peaks at 4.5 GeV/c 2 (see Fig. |J). The latter decays contribute 
a background at values of Mx below the charm mass. They can be identified in part be- 
cause of their lower value of Mxtv compared with fully reconstructed decays. Therefore, 
decays with Mxtv < 3.0 GeV/c 2 were removed. For decays with 3.0 GeV/c 2 < Mxep < 
4.5 GeV/c 2 and Mx < 1.6 GeV/c 2 , the measured hadronic mass Mx was rescaled by 
M B /M X £u- Of the decays in this category that passed the final selections (see Table 1, 
below), this rescaling promoted 46.2% to M x > 1.6 GeV/c 2 in data and 45.8% in simula- 
tion. In order to further remove partially reconstructed D decays, all charged particles in 
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the lepton hemisphere with p > 1.5 GeV/c which were not associated with the secondary 
vertex were tested for their probability to originate at the event primary vertex. Decays 
giving a probability below 0.025 were removed, since this low probability indicates the 
presence of additional secondary particles that were not included. In addition, decays with 
two identified leptons in the same hemisphere were removed, because double semileptonic 
B — > X c £u, X c — > X s iv decays result in a low mass hadronic system and represent a 
background to this analysis. The secondary hadronic system reconstruction and selection 
gave 34583 accepted hemispheres in data and 33769 in simulation. 

Finally, an inclusive search for D* — > Dir was performed. Charged pion candidates 
with 0.4 GeV/c < p < 3.0 GeV/c and p t < 0.7 GeV/c and 7r° candidates with 1.5 GeV 
< E < 3.0 GeV and E t < 0.7 GeV were added in turn to the secondary system X 
and the mass difference AM = Mx-k — Mx was computed. Events with 0.14 GeV/c 2 < 
AM < 0.16 GeV/c 2 and a secondary hadronic mass above 0.6 GeV/c 2 were accepted as 
candidate D* decays. Their mass Mx was then fixed to 2.01 GeV/c 2 . In the simulation, 
these criteria correctly identified 44% of semileptonic decays with a D* meson and Mx < 
1.6 GeV/c 2 . This procedure promoted a further 13.8% to M x > 1-6 GeV/c 2 in data, and 
14.2% in simulation. 



2.4 Boosted Lepton Energy 

In order to improve the separation of B — > X u iv from B — > X c iv and other background 
sources, the lepton energy in the B rest frame was determined. For each decay, the 
energy of the B hadron (see Fig. 3) was estimated as the energy sum of the identified 
lepton, the secondary hadronic system and the neutrino energy. The neutrino energy was 
computed from the missing energy in the hemisphere corrected by a function of the Exi 
energy determined from the simulation [|l7j]. Neutrino energies in the range of 1.5 GeV 



< E v < 25 GeV and a minimum B energy of 25 GeV were required. The resolution of 
the neutrino energy in B — > Xtv decays was estimated to be 3.6 GeV. The resulting 
resolution of the B energy was studied on simulation and found to be 9.8% for 80% of 
all inclusive semileptonic B decays and 15.4% for the remaining decays. The B direction 
was taken as the vector joining the primary vertex to the jet secondary vertex. For those 
decays without a reconstructed secondary vertex, the sum of the momentum vectors of 
the hadronic system, of the lepton and of the missing momentum was computed. The 
angular resolution of the B direction was estimated to be 3.2° for semileptonic B decays. 

The lepton energy E\ was computed in the frame defined by the estimated B energy 
and direction. The resolution on the E\ reconstruction was studied on simulated events 
and found to be 14% for 81% of the selected decays. The values of the resolution obtained 
by analysing B — > X u tv and B — > X c tv events separately were found to be consistent. 



2.5 b — > u enrichment 

A procedure was developed to select separate samples, enriched or depleted in b — > u 
transitions, independently of the reconstructed hadronic mass. It relies on the sign of the 
lepton impact parameter relative to the secondary vertex position and on the presence 
of identified kaons in the same hemisphere as the lepton. For each hemisphere, with a 
reconstructed secondary vertex, the lepton impact parameter d e sec was computed relative 
to this vertex and signed using the lifetime convention, i.e. it was signed negative if the 
lepton appeared to originate between the primary and the secondary vertex, and positive 
if it was downstream of the secondary vertex. In B — > X c lv decays, the secondary 
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vertex corresponds mainly to the charm decay vertex. Consequently, b — > c semileptonic 
transitions tend to give leptons with negatively signed impact parameters, since the lepton 
comes from the B decay vertex. But in B — > X u £v transitions, the secondary vertex 
coincides with the B decay vertex and thus with the lepton production point, so the 
impact parameter signing depends only on resolution effects and is positive or negative 
with equal probability. Therefore, events with a reconstructed secondary vertex and a 
significantly negative lepton impact parameter were assigned to the b — > u depleted class. 
Requiring d l sec < —0.015 cm selected 29.4% of the decays fulfilling the final selection in 
data, and 28.7% in the simulated background sample. Those with a single secondary 
particle not identified as a kaon or a proton were assigned to the b — ► u enriched class. 

The detection of a strange particle in the semileptonic B decay was also used to separate 
cascade B —+ X c £v followed by X c — > K ± X decays from B — > X u £v transitions, where 
the production of strange particles is suppressed because they can originate only from the 
spectator s quark in B s decays or from the production of an ss pair in the hadronisation 
process. In the same hemisphere as the lepton, 36.5% of reconstructed decays contained 
an identified K ± , or proton in the data, and 37.8% in the simulated backgrounds, 
while the simulation predicts 14% in b — > u transitions. These decays were assigned to 
the b — > u depleted class if d sec < 0. 

3 Results 

Candidate semileptonic B decays were further selected by imposing the following selec- 
tion in order to remove background and poorly reconstructed decays. The summed energy 
of the hadronic system and the lepton was required to be larger than 12 GeV and larger 
than 70% of the jet energy. Decays with an invariant mass of the secondary hadronic 
system and of the lepton Mxi below 2.0 GeV/c 2 were also removed. Finally, decays in 
which the lepton charge had a sign equal to that of the hadronic system were discarded. 
These criteria selected 12134 decays in data. The Z° simulated sample, which contained 
no b — > u transitions, gave 11695 expected events, while from a dedicated signal sample 
an efficiency for B — > X u £u, with £ = e, //, of (9.3±0.3)% was obtained. The background 
composition was studied on simulation, for E\ > 0.8 GeV, and found to consist of 90% of 
B — > X c £v decays, 8% of cascade B — > X c — > X s £v decays and 2% of D — > X s £v decays 
and misidentified hadrons. 

Selected decays were divided into four independent classes according to the recon- 
structed secondary hadronic mass Mx and the b — > u enrichment criteria, described in 
section 2.5. These are: i) b — > u enriched decays with Mx < 1.6 GeV/c 2 , ii) b — > u en- 
riched decays with M x > 1.6 GeV/c 2 , iii) b — > u depleted decays with M x < 1.6 GeV/c 2 , 
and iv) b — » u depleted decays with M x > 1.6 GeV/c 2 . The M x value of 1.6 GeV/c 2 
was chosen on the basis of simulation studies. These showed that this value (a) was large 
enough to ensure a reduced model dependence in the extraction of \V u b\ (b) was 

sufficiently below the D mass to suppress the bulk of B — > X c £u decays, and (c) minimised 
the statistical error on |V^&|/|Vd,|. 

The numbers of events selected in data, the numbers of expected background events, 
and the expected fractions of the total number of signal events in the four classes are 
summarised in Table |TJ. The background was rescaled by the normalisation factor obtained 
from the fit described in section 3.1 including the ±0.01 error on this normalisation. For 
decays selected in the low Mx and b — > u enriched class, which is expected to contain 
68% of the b — > u signal, an excess of 214±56 events above the expected background was 
found in the data (see Figs. No significant excess was observed in the other classes, 
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where the y 2 probability of a deviation of the data from the prediction larger than that 
observed is 30%. 

As a cross-check, the analysis was repeated using both anti-6-tagged events and decays 
with same-sign lepton and hadronic vertex combinations. All other selection criteria were 
kept as in the main analysis. Both these samples are expected to be depleted in signal 
b —> u decays but they are sensitive to possible discrepancies between data and simulation 
in the description of backgrounds. The numbers of selected decays in the low Mx and 
b — > u class were computed. For the anti 6-tagged sample 32 events were observed in the 
data compared to 33±3 expected from background. The same-sign sample consisted of 
340 events in data with 317±8 expected from backgrounds. No excess of events in data 
was observed in either of these samples. 

Table 1: Numbers of events selected in the data, expected background events, and ex- 
pected fractions of the total number of signal events. The background is rescaled by the 
normalisation factor obtained from the fit described in section 3.1, including its ±0.01 
error 



Selection 


b — > u enriched 


b — > u depleted 




Data 


Back. 


Sig. 


Data 


Back. Sig. 


M x < 1.6 GeV/c 2 


2292 


2078±30 


68% 


1081 


1118±19 


9% 


M x > 1.6 GeV/c 2 


5017 


5019±60 


17% 


3744 


3618±46 


6% 



In order to check the interpretation of the excess of events in the low Mx and b — > u 
enriched sample as a signal for B — > X u lv transitions, a search for decays into the B — > tt£u 
and B — > ptv exclusive final states was performed. Decays with 1.0 GeV < E\ < 3.0 GeV 
and a reconstructed hadronic system consisting either of a single charged particle or of 
two particles with total charge Q = or ±1 were selected. The lepton energy requirement 
further suppressed the non B — > X c £u backgrounds. The Mx distribution shows an excess 
of events in the data compared to the expected backgrounds in good agreement with the 
expectation from B — > i\iv and B — ► piv exclusive final states (Fig. 8). As a cross-check, 
the analysis was repeated for same-sign combinations of the lepton and the hadronic 
system. This class receives signal contributions only from partially reconstructed decays 
like B + — > 7r + (7r _ )£ + z/, where the (tt~) is not reconstructed, and background events. No 
significant excess of events was observed (see Fig. 8). 



3.1 Extraction of |VU|/|K*| 

The numbers of events in each decay class and their E\ distributions were used to deter- 
mine the value of |Kij,|/|V^,| by a simultaneous binned maximum-likelihood fit. The ratio 
|Kife|/|Kb| is given by the ratio of X u tv to X c tv decays through the relationship [I8] , |I9"|] : 



[Khj _ 0.00445 ( BR(B^X u ep e ) O105\ 1/2 

|V ci ,| 0.04110 A \BR{B^X c e& e ) A 0.0O2J (1) 

x(l±0.055 QCD ±0.015 m J 

In the fit, the overall data to simulation normalisation and the value of |V^6|/|Vrf,| were 
left free to vary while the non B — > X c lv backgrounds were kept fixed to the fractions 
predicted from the simulation. This significantly reduced the systematic uncertainties 
from the lepton identification and other sources by absorbing their effects in this overall 
normalisation factor. 
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The result of the fit (see Fig. 7) was: 

\v ub \l\v cb \ = 0.103 ±8:8S (stat.) 

with the normalisation factor 1.013±0.011. 

3.2 Stability checks 

The stability of the fitted value of |T4b|/|Vd,| was checked in various ways. Repeat- 
ing the analysis for electrons and muons separately gave 0.095±0.017 and 0.107±0.014 
respectively. Moving the cut on M x from 1.6 GeV/c 2 to 1.05 GeV/c 2 or 1.90 GeV/c 2 , 
which changed the signal-to-background ratio from 0.10 to 0.15 or 0.055, changed the fit 
result to IKj&I/I^c&I = 0.098±0.013 or 0.105±0.013, respectively. Several other selection 
criteria were also varied or dropped, including the cuts on event ^-tagging, lepton p t , E X t, 
Mxe, Eb and the secondary vertex decay distance significance, and the results were found 
to agree within the errors. Removing the scaling of the hadronic invariant mass Mx, de- 
scribed in section 2.3, gave 0.105±0.013. Excluding the information on the lepton energy 
gave 0.106±0.015. Keeping the simulation to data normalisation fixed gave |V^|/|\4b| = 
0.108±0.009. 

Finally, the analysis was repeated with an improved rejection of double semileptonic 
decays B — > X c £u, X c — > X s £v. Because of the presence of an additional neutrino 
and the lower secondary charged particle multiplicity, these decays can be confused with 
B — > X u lv transitions. A loose e and \x identification procedure, based on the dE/dx 
measured in the TPC, the measured ratio of calorimetric energy to momentum, E/p, 
and the hit pattern in the Hadron Calorimeter, was developed. This had an efficiency of 
75% for leptons not already identified by the standard tagging procedure and a hadron 
misidentification probability of 18% in the momentum range 3 GeV/c < p < 10 GeV/c. 
Candidate decays were rejected if the secondary charged particle multiplicity (not includ- 
ing the seed high p t lepton) did not exceed two and at least one secondary particle was 
tagged as a lepton with these looser requirements. This gave iK^I/jV^I = 0.100±0.013. 

3.3 Systematic Uncertainties 

Four categories of systematic uncertainties were considered. The first three affect 
primarily the estimation of the large b — *> c backgrounds. The fourth affects primarily 
the evaluation of |V^&|/|V^,| from the observed excess in the b — > w-enriched sample. The 
results for IV^I/IK&I are summarised in Table ||. 

3.3.1 Systematic errors from charm decays 

The description of charm decays affects the fraction of B — > X c lv transitions that 
were accepted in the b — ► u enriched sample. First, the branching fractions for D decays 
into final states with low charged multiplicity were considered. These are D — ► K°X, 
which in the simulation contributes 40% of the background from b —>■ c decays in the 
b —>■ u enriched sample, and D decays in and 1 prong final states, which contribute 46%. 
Varying their branching ratios within the uncertainties of the Mark III measurement [pOfl 
changed the fitted value of |Vub|/|K&| by ±0.0062 and ±0.0025, respectively. Varying the 
ratio of prompt B semileptonic decays to cascade and charm decays by the uncertainty 
on the ratio BR(6 -> £)/BR(c -> i) = (0.106±0.002)/(0.098±0.003) J2J contributed an 
uncertainty of ±0.0042. 
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3.3.2 Uncertainties in B hadron production and decay 

The first component is due to the fraction of B s and beauty baryons produced. Due 
to the rejection of kaons and protons associated with the lepton hemisphere these beauty 
hadrons do not significantly contribute to the b — > u enriched sample. Propagating 
the uncertainty of ± 0.020 on the sum of B u and B^ meson fractions ]19] contributed 



0.0039. Varying the parameter in the Peterson b fragmentation function according to 
the uncertainty in the fraction of the beam energy taken by the beauty hadron, < Xb > 
= 0.702±0.008, contributed 0.0010. Varying the inclusive b lifetime by the uncertainty 
of the present world average, (1.564±0.014) ps contributed 0.0011. The branching 
fraction for inclusive double charm production in B decays was fixed to 0.15 and varied by 
±0.03. This changed the fit result by ±0.0025. Finally, the dependence on the production 
rate of D* and D** mesons (where D** denotes either a non-resonant D*n final state 
or a Z)j*' ) higher excited charmed meson state) in B semileptonic decays was studied. 
These states flip the charge of the resulting D meson and thus increase the charged 
multiplicity in b hadronic decays. Using the values BR(B -> D*£u) = (0.046±0.003) @ 
and BR{B -> D**£u) = (0.034±0.006) @, gave a sum of 0.080±0.007, corresponding to 
a contribution to the systematic error of 0.0033. In addition, the amount of non-resonant 
7r states was varied from zero to 50% of the D** yield. This contributes systematic 
uncertainties both from the shape of the lepton energy spectrum and from the vertex 
topology and secondary charged multiplicity corresponding to an uncertainty of ±0.0065. 

Finally the systematic error from the model of the shape of the lepton spectrum in the 
background b —> c semi-leptonic transitions was estimated. The lepton spectra observed 
for the three signal-depleted classes were found to agree with those from the simulation 
(see Fig. 6). In addition, a sample enriched in B — ► D*(X)£v decays was compared with 
the simulation prediction and was also found to be in good agreement. The spectrum 
predicted by the Delphi simulation program was compared with that from the ISGW- 
2 model implemented in the EvtGen decay generator [[25j . The relative contributions 
of the different charm states were set to the central values discussed above in order to 
be sensitive only to the difference in the predicted shape of the lepton spectrum. The 
difference of 0.0020 in the fit result was taken as the corresponding contribution to the 
systematic uncertainty. 



3.3.3 Detector dependent systematics 

The first source is due to the efficiency and purity of the lepton identification. Efficiency 
and misidentification probability for muon and electron tagging were extracted from both 
the simulation and data as discussed in 2.2. The central values for the simulation were 
changed within these errors and the corresponding changes on |Kib|/|K*| were found to 
be 0.0015 and 0.0020, respectively. 

The second component of detector systematics is due to the hadronic mass and rest- 
frame lepton energy resolution. The former depends on both the secondary hadronic 
multiplicity and the single particle energy resolution. The secondary hadronic multiplicity 
was studied on the b — > u depleted sample to avoid possible biases from the presence of 
signal events in the data. The total (charged) multiplicity of the secondary hadronic 
system was measured to be 3.64±0.01 (2.74±0.01) in the data and 3.63±0.01 (2.72±0.01) 
in the simulation. The systematic uncertainty was evaluated by considering a change of 
±0.02 units of multiplicity. This gave a variation of | V c b\ of ±0.0065. The component 
due to the particle energy resolution is dominated by the resolution of neutral particles. 
The systematic effects were checked by decreasing by 2% the resolution on Mx for those 
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decays with neutral particles. This range was defined by a x 2 analysis of the hadronic 
mass in real and simulated b — > u depleted decays. It also corresponds to the effect from 
the electromagnetic energy resolution typically measured in Z° decays. The corresponding 
systematic uncertainty was found to be 0.0010. The resolution on the neutrino energy 
reconstruction was varied by 10% and the effect was propagated to the E\ resolution, 
giving a systematic error contribution of 0.0010. 

Thirdly, the possible systematics in the decay classification were studied. These system- 
atics depend both on the measurement and sign of the lepton impact parameter relative 
to the secondary vertex and on the kaon identification. The lepton impact parameter 
systematic has two components. The first is due to the lepton extrapolation and the 
second to the secondary vertex reconstruction. The effect of changing the resolution on 
the lepton track extrapolation was computed by smearing the resolution on the lepton 
impact parameter by 5%, which corresponds to the maximum discrepancy observed in 
the resolution functions obtained in data and in simulation ||26|| . The component due 



to the secondary vertex position was evaluated by smearing the resolution on its decay 
length by 50 /im, which corresponds to the additional smearing that increases by one 
the x 2 °f a data-to-simulation comparison of the decay distance distributions. Summing 
these two effects in quadrature gave a systematic uncertainty of ±0.0070. The kaon tag- 
ging efficiency was varied by ±2.5%, corresponding to the largest observed deviation of 
the performance of the hadron identification tagging in data and simulation, and the 
corresponding uncertainty on the result of the fit was 0.0025. 

Finally, the statistical error on the efficiency for selecting signal events contributes 
±0.0015 to the systematic uncertainty. 



3.3.4 Uncertainties in the B — > X u £u model 



The predicted shape of the invariant-mass distribution in the B — > X u iv decay, de- 
pends mainly on the kinematics of the heavy and spectator quarks inside the B hadron 
and on the b quark mass. Further, the hadronisation process, transforming the uq system 
into the observable hadronic final state, represents an additional source of model uncer- 
tainties. These uncertainties were studied using a dedicated generator that implements 
different prescriptions for the initial state kinematics and the resonance decomposition of 
the hadronic final states ||13|| . 

Varying the b quark mass by ±100 MeV/c 2 [18|,p7,p8|, gave a systematic error of 
±0.0045. The value of the b quark mass introduces also an uncertainty in the extraction 
of |14&|/|Vy from the observed B — ► X u tv rate |T8|j27[| , see Eq. 1. This gives a total error 
contribution of 0.0047. The average kinetic energy of the b quark in the hadron, < p\ >, 
has been evaluated both from theory and from fits to experimental data. Results are 
scheme or model dependent and depend on the method used in their derivation, but point 
to the value of the parameter /x 2 = (0.5±0.1) GeV 2 . This variation contributes ±0.0015 
from the uncertainty on the hadronic mass spectrum and ±0.0024 from the derivation of 
the ratio | V u b\/\ V c b\, giving a total of ±0.0028. The description of the motion of the b quark 
inside the heavy hadron also contributes uncertainties. The momentum distributions from 
the ACCMM model P^ ], a shape function PD|jn| parametrised as f(z) = z a {l — cz)e~ cz , 
and the parton model [^] were compared. For the ACCM, the b quark pole mass was 
kept fixed at its central value of mj = 4.82 GeV/c 2 and the pp value and the a and c 
coefficients in the QCD structure function were chosen to reproduce the same value of 
< p\ >. For the parton model, the Peterson form |33| of the fragmentation function was 
adopted with e& = 0.0040. A systematic uncertainty of ±0.0025 was evaluated. 
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Table 2: Summary of the systematic uncertainties on \V u b\/\V c b 



Source 


Value ± Range 


Syst. Error 


BR(D -> K^X) 


0.53 ±0.05 


0.0062 


~BR(D -> 0, 1 prong) 


0.22 ± 0.02 


0.0025 


B(6 -> £)/BR(c -> £) 


0.106±0.002 
0.098±0.003 


0.0042 


Charm decay sub-total 




0.0079 


f(B u )+f{B d ) 


0.802 ±0.020 


0.0039 


<x b > 


0.702 ±0.008 


0.0010 


b lifetime 


1.564 ±0.014 


0.0011 


BR(6 -> ccs) 


0.15 ±0.03 


0.0025 


BR(5 -> L>*£z/ + 


0.080 ±0.007 


0.0033 


BR(5 -> D(*)vr£z//BR( J B -> 


0.25 ±0.25 


0.0065 


Lepton spectrum shape 




0.0020 


B production and decay sub-total 




0.0090 


e / /i id. Efficiency 


±2.5% 


0.0015 


e / fi id. Purity 


±10% / ±4.5% 


0.0020 


Hadronic Multiplicity 




0.0065 


Neutral Energy Resolution 




0.0010 


Missing Energy Resolution 




0.0010 


t Impact Parameter 




0.0070 


K id. Efficiency 


±2.5% 


0.0025 


Signal Efficiency 


± 3.0% 


0.0015 


Detector-dependent sub-total 




0.0104 


m b 


4.82 ±0.10 


0.0047 


<p 2 b > 


0.5 ±0.1 


0.0028 


b Kinematic Model 




0.0025 


Hadronisation Model 




0.0060 


QCD corrections 




0.0050 


b —>■ u model sub-total 




0.0099 


Total 




0.0186 



The production of the hadronic final states from the uq pair was simulated according 
to both the ISGW-2 [24| exclusive and a fully inclusive model based on parton shower 



fragmentation |Tl|. The ISGW-2 model approximates the inclusive B — ► X u lv decay 
width by the sum over resonant final states, taking into account leading corrections to the 
heavy quark symmetry limit. The predicted branching ratios for the different resonant 
final states were used to define the hadronic system emitted with the lepton. Another 
approach is to assume that, at sufficiently large recoil u quark energies, the uq system 
moves away fast enough to resemble the evolution of a jet initiated by a light quark q in 
e + e~ — > qq annihilation. This was simulated by first arranging the uq system in a string 
configuration and then requiring it to fragment according to the parton shower model. 
Due to the extreme assumptions of the two models adopted, the resulting difference of 
0.0160 in the fitted value of was assumed to correspond to a 90% confidence 

region and the ±la syst was estimated to be ±0.0060. 



11 



Additional sources of theoretical systematics, arising from the perturbative part of the 
evaluation and the contribution of non perturbative corrections of order l/m\ contribute 
a ±0.0050 systematic error 



4 Summary and Discussion 



The value of the ratio |T46|/|K*| was measured using a novel technique. The technique 
uses the reconstructed mass Mx of the secondary hadronic system produced in association 
with an identified lepton in the semi-leptonic decay of a B hadron and the rest-frame 
energy spectrum of that lepton. The b —>■ u signal is enriched using identified kaons and 
protons and the lepton impact parameter with respect to the secondary vertex. The result 
obtained is 

\V ub \/\V cb \ = 0.103 toill (stat.)±0.016 (syst.) ± 0.010 (model). 

Here the systematic error quoted combines the charm decay, B production and decay, and 
detector-dependent systematics evaluated above, which all primarily affect the estimation 
of the b — > c backgrounds, and the b — > u model uncertainty is quoted separately. The 
technique adopted in this analysis allowed both sources of systematic uncertainty to be 
reduced. At the critical points of the analysis, the behaviour of the data agrees well with 
the expectations from the simulation. The result is found to be stable with respect to 
variations in the analysis procedure. 

There remains a possible further model dependence arising from a biased sampling 
of the decay phase space in B — > X u lv transitions. The stability of the result when 
the M x cut was moved from 1.6 GeV/c 2 down to 1.05 GeV/c 2 or up to 1.90 GeV/c 2 
(section 3.2) argues against this. To further check this possibility, the relative weights 
of different regions in the Mx — El plane were analysed. Four regions of the Mx — E\ 
plane were defined by selecting decays with Mx and E\ above and below 0.8 GeV/c 2 
and 1.75 GeV, respectively. The fit was repeated separately for these four regions. The 
statistical weights of the four regions are given in Table [3]. The result agrees with the 
expectation of a higher contribution from charmless semileptonic B decays in the low- 
mass, high-energy and high-mass, low-energy regions and indicates no strong bias in the 
weighting of the decay phase space. 



Table 3: The weights of the four regions in the Mx-E\ plane, used to check the phase 
space sampling in the determination of \ V u b\/\V c b\. 





0.1< M x <0.8 GeV/c 2 


0.8< M x <1.6 GeV/c 2 


1.75< El <3.0 GeV 


0.33 


0.21 


0.1< E\ <1.75 GeV 


0.15 


0.29 



While this analysis extracted the ratio of CKM elements from the fitted 

fraction of candidate B —>■ X u lv decays, it is also interesting to extract explicitly the 
the charmless semileptonic branching ratio. This was obtained from the fitted result, 
assuming \V ch \ = (38.4±3.3) xl0~ 3 and r b = (1.564±0.014) ps |H§. The result was: 

BR(B -> XJu) = (1.57±0.35 (stat.)±0.48 (syst.)±0.20 (\V cb \) ± 0.01 (r 6 ) 
± 0.27 (model)) xl0~ 3 x ( | Kfol /0.0384) 2 . 

where the contribution of correlated model systematics in the derivation of \V ub \ and \ V cb \ 
was taken into account. 
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Figure 1: The invariant mass spectrum Mx of the u-spectator quark system for inclusive 
B — > X u tv decays at the parton level obtained with a dedicated decay generator jjT3J 
(upper plot) and the fraction of the decays with Mx below a given value (lower plot). 
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Figure 2: The total (upper plot ) and charged (lower plot ) multiplicity of the reconstructed 
secondary hadronic system in selected decays for data (points with error bars ) and simu- 
lation (histogram). 



16 




10 20 30 40 10 20 30 40 



Missing Energy (GeV) X 1 v Energy (GeV) 



Figure 3: Energy Distributions: i) energy of the hadronic system (upper left plot), ii) 
energy of the hadronic system plus the lepton (upper right plot), in) missing energy 
(lower left plot) and iv) reconstructed B energy (lower right plot). The dots with error 
bars represent the data and the histograms the simulation. 




Figure 4: Invariant mass Mxiv of the reconstructed B decay: data (points with error 
bars) and simulated background (histogram) (upper plot) and b — > u simulated signal 
(lower plot). 
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Figure 5: Invariant mass Mx of the reconstructed secondary hadronic system in selected 
decays for data (points with error bars) and simulation (histogram). The plots show the 
b — > u depleted sample (upper plot), the b — ► u enriched sample (medium plot) and the 
b — > u signal (lower plot). The vertical lines correspond to the value chosen for the low 
M x selection. 
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Figure 6: The E\ distribution for the decays in the four selected classes: i) b — > u 
enriched decays with M x < 1.6 GeV/c 2 (upper plot), ii) b — ► u enriched decays with 
M x > 1.6 GeV/c 2 (lower left), in) b — > u depleted decays with M x < 1.6 GeV/c 2 (lower 
central plot), and iv) b — > u depleted decays with M x > 1.6 GeV/c 2 (lower right plot). 
Data are indicated by the points with error bars, the b — ► X u lv signal by the dark shaded 
histograms, the b — > X c iv background by the medium shaded histograms, and the other 
backgrounds by the light shaded histograms. 
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Figure 7: Background subtracted El distributions: the b — » u enriched decays with Mx < 
1.6 GeV /c 2 (upper plot) and b — > u depleted decays with M x < 1.6 GeV/c 2 (lower plot). 
The background was rescaled by the fitted normalisation factor. The shaded histograms 
show the expected E\ distribution for signal B — ► X u tv decays normalised to the amount 
of signal corresponding to the fitted \ V u b\/\V c b\ value. 
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Figure 8: Invariant mass Mx distributions for b — > u enriched decays with a secondary 
hadronic system consisting of either a single charged particle or two particles forming a 
neutral or unit-charge secondary system for opposite sign and neutral- charge (left plots) 
and same-sign (right plots) hadronic-lepton system. The upper histograms show the 
expected distribution from backgrounds, the points with error bars the data. In the lower 
histograms the background subtracted data are compared with the expected distribution 
from signal B — > X u lv events. 



